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Abstract. The concept of time reversal (TR) of scalar wave is reexamined from basic 
principles. 

Five different time reversal mirrors (TRM) are introduced and their relations are ana- 
lyzed. For the boundary behavior, it is shown that for paraxial wave only the monopole 
TR scheme satisfies the exact boundary condition while for the spherical wave only the 
MD-mode TR scheme satisfies the exact boundary condition. 

The asymptotic analysis of the near-field focusing property is presented for two dimen- 
sions and three dimensions. It is shown that to have a subwavelength focal spot the TRM 
should consist of dipole transducers. The transverse resolution of the dipole TRM is linearly 
proportional to the distance between the point source and TRM. The mixed mode TRM 
has the similar (linear) behavior in three dimensions but in two dimensions the transverse 
resolution behaves as the square-root of the distance between the point source and TRM. 
The monopole TRM is ineffective to focus below wavelength. 

Contrary to the matched field processing and the phase processor, both of which resemble 
TR, TR in a weak- or non-scattering medium is usually biased in the longitudinal direction, 
especially when TR is carried out on a single plane with a finite aperture. This is true for 
all five TR schemes. On the other hand, the TR focal spot has been shown repeatedly in the 
literature, both theoretically and experimentally, to be centered at the source point when 
the medium is multiply scattering. A reconciliation of the two seemingly conflicting results 
is found in the random fluctuations in the intensity of the Green function for a multiply 
scattering medium and the notion of scattering-enlarged effective aperture. 



Time reversal (TR) is the process of recording the signal, time- reversing and re-propagating 
the signal. When the signal is from a localized source the time reversed field is expected to 
focus on the source. Time reversal of acoustic waves has led to applications in ultrasound 
and underwater acoustics including brain therapy, lithotripsy, nondestructive testing and 
telecommunications [17] . An even greater potential holds for the time reversal of electro- 
magnetic waves which is closely related to optical phase conjugation [28J. 




Recently, time reversal experiments with optical and micro waves demonstrate phase con- 
jugation of optical near field and the robustness of turbidity suppression by optical phase 



Motivated by these exciting experiments, we reconsider time reversal from basic principle, 
discuss various TR schemes and their mutual relations, and make general observations about 
the focusing properties of near-field and far-field time reversal. 

First, we review the principle of perfect time reversal in a closed cavity based on the Green 
second identity and discuss the focusing properties of the perfect TR kernel (the Porter- 
Bojarski kernel). We review the well known fact that in the free space the TR focal spot of 
a field with source is always diffraction-limited even when all the evanescent components are 
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recorded by time reversal mirror (Section [2] and |4j). The perfect time reversal is "perfect" 
only for the source-free fields. On the contrary, any medium inhomogeneities, no matter how 
weak, always give rise to a nonvanishing evanescent component in the time reversed field. 

Next we derive four other time reversal schemes (monopole, dipole and two mixed modes) 
from the Porter-Bo jarski kernel and analyze their behaviors with one planar TRM. We derive 
their mutual relationships for the paraxial wave. For the boundary behavior, we show that 
for paraxial wave only the monopole TR scheme satisfies the exact boundary condition while 
for the spherical wave only the MD-mode TR scheme satisfies the exact boundary condition. 

A main interest here is the near-field focusing property of various TRMs. We show that 
to achieve a subwavelength focusing the near-field time reversal mirror should involve dipole 
fields in the sense that the TRM records not just the field but the normal gradient of the 
field and/or back-propagate the dipole field in proportion to the phase-conjugate recorded 
data. The dipole TRM records and re-transmits the dipole field while the mixed mode TRM 
involves only one of the processes. The monopole TRM records and transmits only the 
monopole field. The transverse resolution of the dipole TRM is linearly proportional to the 
distance between the point source and TRM. The mixed mode TRM has the similar (linear) 
behavior in three dimensions but in two dimensions the transverse resolution behaves as 
the square-root of the distance between the point source and TRM. The monopole TRM is 
ineffective to focus below wavelength (Section [3]). 

We also point out that TR focal spot with a finite aperture is generally closer to the TRM 
than the source point in the free space or a weakly scattering medium such as consisting 
of phase scatterers (Section [5]). In other words, in the absence of multiple scattering in 
the medium, TR as focusing or imaging method is biased. In comparison, the conventional 
matched field processor and the phase processor always produce a centered focal spot (Section 
5]). Lastly in Section |6j we explain how the random effect in multiple scattering can restore 
the centeredness of TR focal spot on the source point as well as the stability condition for 
the robustness of turbidity suppression by phase conjugation [27] . 

2. Perfect time reversal 
Consider a real-valed signal in the temporal Fourier representation 



U(r,t) = J u(r,uj)e~ iuJt du. 



The real-valuedness of U implies that u(t,—u) = u*(t,u). Therefore, time reversal is 
equivalent to phase conjugation of every frequency component u(r,u) — > u*(r,u), Vcj, r. 
Ideally, TR turns a divergent wave into a convergent wave, thus focusing wave energy. 

A practical way of phase conjugation within some control domain is to operate TR from 
the boundary of the control domain as we consider next. 

Consider a monochromatic scalar wave u propagating in a medium characterized by the 
refractive index n which in general varies in space. Such a u satisfies the Helmholtz equation 

(1) Au + k 2 n 2 u = f 

subject to suitable boundary conditions where / is the source of compact, localized support. 
We assume that the medium is dissipationless, i.e. n 2 is real-valued. 

In time reversal, the source / emits a wave field u which is then recorded at the boundary 
dD of the (topologically open) domain D. The typical domain can be one of the following 
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Figure 1. Three representative domains: (i) a bounded domain (left), (ii) an 
infinite slab (middle) and (iii) a half space. The dark areas represent medium 
inhomogeneities (scatterers). 



three kinds: bounded with a closed boundary, infinite slab with two planar boundaries, or 
half space with one planar boundary (Figure [l]). 

Consider a bounded domain D first. Let both u and du/dn be recorded at dD, phase- 
conjugated and back-propagated into the domain D and let v(r) be the resulting wave field. 
Let G(x, y) be the Green function of 0. The time reversal principle originally proposed in 
[H] is described by 



(2) 



vir) 



a 



dn 



-der(r') 



dD 



du*{r') 
dn 



G(r,T')da(r') 



Clearly, v(r) is a source-free solution of the Helmholtz equation in D. 
Whenever Sommerfeld's radiation condition 



(3) 



d_ 
Or 



G(r,r') w ikG(r,r') 



for all r' sufficiently far away from r, ^ can be easily extended to the second type of 
domains (infinite slabs) through a limiting procedure. Under this condition there is no 
difference between type (i) (bounded) and type (ii) (infinite slab) domains. 

Let A be the time reversal mirror (TRM) considered as a part of dD. At each point 
on the boundary dD\A we impose either the sound soft boundary condition (G = 0) or 
the sound hard boundary condition (dG/dn = 0). The boundary condition on TRM A is 
flexible except that the presence of TRM does not affect the near-field behavior of the Green 
function, see (26)-(27). We call either type (i) or (ii) domain of Figure [l] closed cavity. 



By the second Green's identity we have 



*/ r\dG(r, r') . 

u (r ) da(r ) 

dD on 



du*(r' 



3D 



dn 



■G(r,r')da(r') 



(4) 



/ u*(r')(A + k 2 n 2 (r'))G(r,r')dr' - [ (A + k 2 n 2 (r'))u* (r')G(r , r')dr' . 

J D J D 



The time reversal operation described by ^ is perfect for the field inside D when the field 
is source- free inside D (i.e. u satisfies Q with f\ D = and the second integral on the right 
hand side of Q drops out) 

(5) v(r) = u*(t), reD 

while v(r) = 0,Vr ^ D by Q. Hence (@) has been proposed as the basis for time reversal in 
a closed cavity [HI E] • 

On the other hand, if D contains sources as in typical applications, then the time reversal 
described by Q is diffraction-limited. Indeed, 



(6) v(r) = u*(r)- f*(r')G{r,r')dr' 

Jd 

= [ [G*(v,r')-G(v,r')]r(v')dr' 
Jd 

= [ K FB (r,r')f*(r')dr', reD 
Jd 

with the Porter-Bojarski kernel 

(7) K PB (r,r / ) = -2S[G(r,r / )] 

where 3= denotes the imaginary part, c.f. [8]. With a real-valued n 2 , Kpv(r, r') is an example 
of a source- free radiation field in r for each r'. In the free space, ([7]) is always diffraction- 
limited (see Section [4]). 

The reason that Kpb 7^ G* is because the time symmetry is broken when the field is phase 
conjugated but not the source. A time-reversed source is a sink. Suppose now the point 
source becomes a point sink after emitting the initial wave and suppose the sink can absorb 
the time-reversed, incoming wave and prevent the re-emission of the outgoing wave with 
efficiency p > 0, [ID]. Then @ becomes 

(8) v(r) = G*(r,r')-(l-p)G(r,r ), reD, 

which, denoted by K$, is the TR kernel with a point sink of efficiency p. When p = 1 
then Ks = G*. As we shall see in Section [4], the TR focusing property for Kg is infinitely 
enhanced over that for 

3. MONOPOLE AND DIPOLE TRMS 

In this section we derive approximate TR schemes from ^ using Sommerfeld's radiation 
condition 

For a closed cavity, using the sound soft/hard condition on dD\A, we write 



JdD on J dD dn 

(r )dr 

which yields the alternative expression for the Porter-Bojarski kernel for a closed cavity 
(9) K PB (r,r") = J G\v\v") d -^^da(v>) - J G(v,v')da(v') 



(r ' r ^ dn (l " ^ ~ / dn (r ' r ) (r ^ 
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Now we pretend |r' — r"|, |r — r'| 3> 1 and use condition ([3| to motivate two other TR 
schemes. We consider all three types of domains depicted in Figure [TJ 

Replacing dG/dr with ikG or ikG with dG/dr in ^ leads to, respectively, the monopole 
TR kernel 

(10) K M {r,r") = -2ik [ G*{r',r")G{r,r')da{r') 

J A 

and the dipole TR kernel 

(11) K D (r,r") = -j J |-G*(r',r")|^G(r,r')^(r'). 



(12) 
(13) 



In addition, we introduce mixed mode TR kernels 

, ,„dG(v,r') 



K M D(r,r") 
K DM (r,r") 



da(r') 



dG^yy) 

A ^—G(r,r)da(r) 



which are parts of 

Analogously we define the various TR schemes for a source- free field u: 



(14) 


VM 


(r) 


= -2ik [ u*{r')G(r,r')da(r') 

J A 


(15) 




(r) 




(16) 


vmd 


(r) 




(17) 


vdm 


(r) 


-- Z^(-w) 


(18) 


VPB 


(r) 


= %D(r) + ^DM(r) 


The various 


TR schemes 


can 


be readily implemented for the acoustic 



TRM, the transducers record the normal pressure gradient and emit the normal dipole field in 
proportion to the phase-conjugate recorded data. For the mixed mode TRM, the transducers 
record the pressure (resp. pressure gradient) and emit the normal dipole (resp. monopole) 
field in proportion to the phase-conjugate recorded data (see, for example, [12] for a similar 
proposal). 

The time reversal schemes described by (|9) 



-(13) and (14h-(18) are main objects of sub- 



sequent analysis. Note that the expressions (9|-(13), ( 14 )-( 18 ) are source-free fields in the 
domain D regardless whether the initial wave is source-free or not. 



3.1. Paraxial wave. In the case of paraxial wave, @, ([10]), g2) and Q are explicitly 
connected with one another as follows. 

Consider the half space D = {z > 0} with a planar TRM on the transverse plane z = 0. 
Let G(z, x, z', x') be the paraxial Green function where z,z' are the longitudinal coordinate 
and x, x' are the transverse coordinates (see Appendix A). Let r = (z , x ) be the location of 
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the source. Assume that the TRM is away from the medium inhomogeneities i.e. n(0, x) = 1. 
Then it can be shown that 

1 i f 

(19^md(^,x,2;o,xo) = -K M (z, x, z Q , x ) + — / V T G*(z ,x Q , 0,x') • V T G(z,x, 0,x')dx' 

= K DM (z,x,z ,x ) 
where denotes the transverse gradient w.r.t. x' (Appendix A). Therefore, 

(20) K MD (r,r ) = K DM (r,r ) = ^K PB (r,r ). 
Likewise, we have 

1 i f 

%d(^,x) = -%(z,x) + — / V t m*(0,x') • V T G(z,x,0,x')cix' 

3.2. Boundary value. Though not fundamental to many time reversal applications, it is 
sometimes desirable for the TR schemes to satisfy the boundary condition 

(21) K x (r,r ) = G'*(r,r ), v x (t)=u*(t), reA 

where X = M, D, MD, DM, PB. We have already seen that in a closed cavity vpb satisfies 
the boundary condition but K PB does not. 

In Appendix [B] we show that, surprisingly, for the paraxial wave, only the monopole TR 
scheme satisfies the boundary condition (21) and for the spherical wave in the free space, 
only the mixed mode TR scheme %d,-^md, multiplied by 2, satisfy (21). 

4. Near-Field focusing 

In this section, we analyze the near-field focusing property of various TR schemes. We 
show that Kd , -Kmd , -Kdm gi ye rise to a subwavelength focal spot of size proportional to 
the distance from the source to TRM while Kyi and Kp B can hardly achieve subwavelength 
focusing. In the free space, the focal spot size for the paraxial wave is always comparable to 
the Fresnel length \f\z$ which is not applicable for zq A. 

Consider the Porter-Bo jarski kernel for the free space Green function in three dimensions, 

is sin (fc|r — r 7 |) 

Z7i\r — r | 

and compare it with the real part of K$ defined by pj) 



(23) m,] = p ""™ - ? ' ">° 

47r|r — r | 

which dominates over the imaginary part for p w 1 and |r — r'| < A/4. 

Resolution of an imaging system generally refers to either the focal spot size or the min- 
imum resolvable distance between two points or lines (two-point or -line resolution). And 
these two ideas are related: the focal spot size is a crude estimate of the minimum resolvable 
distance. 

With all the insights into superresolution techniques, it is now well accepted that in the 
absence of noise and with perfect knowledge of the imaging system there is no limitation 
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to the two-point resolution. Fundamentally, this is because an image with infinite signal-to- 
noise ratio can convey an arbitrary amount of information [3|. 

Since we do not account for the noise explicitly we will focus on the idea of focal spot 
size as a measure of focusing and image sharpness. According to Rayleigh's criterion which 
essentially defines the spot size as the distance of the first zeros to the maximum point, (22) 



has a spot size of A/2 and (23) has a spot size of A/4. But in reality the focal spot is in some 



sense "in the eye of the viewer" and one readily recognizes the vastly sharper graph of ( 23 ) 



than that of (22) 



This contrast is better described by the standard (Houston) criterion used in astronomy 
which is to adopt the "full width at half maximum" (FWHM) [T5]. The FWHM of d22l is 



about A/2 while the FWHM of (23) is zero since the maximum of (23) is infinity (therefore 



infinitely sharp). Indeed, the singularity is the signature of a point source. 

The subwavelength, superresolution of (23) can be understood from Weyl's angular spec- 
trum of plane waves. For the imaging direction z, the free space Green function in three 
dimensions has the angular spectrum representation 



(24) 
where 



ik 
27 



dsids^ 

S3 



exp [ik{si{x - x ) + s 2 {y - y ) + s 3 \z - z \)) 



y/(x- Xq) 2 + (y- y ) 2 + (z- z o y 



s 3 



S3 



1 _ s 2 _ 2 
1 &i 5 2 , 



s\ + s 2 2 < 1 



s? + si - 1, si + si > l 



The integrand in (24) with sf + s 2 , < 1 corresponds to the propagating (homogeneous) 



wave and that with s'j + s 2 , > 1 corresponds to the evanescent (inhomogeneous) wave. Let Wb 
and We denote the homogeneous and the evanescent fields, respectively. Since k\J s\ + s| is 
the transverse wavenumber, Wn contains the diffraction-limited information about the point 
source while We contains the sub-wavelength information about the point source. 

The first important observation is that We is always real-valued. Consequently, the TR 
kernel ifpe in the case of G = Gq contains no evanescent component. Secondly, Wn is 
bounded and in the 5-neighborhood of the source point the estimate holds: |9£[Wh]| = 0(5). 
In other words, in the immediate vicinity of the source point 3? [Go] essentially coincides with 
We- This means that 9ft [Go] contains all the subwavelength information about the point 
source. 

Therefore, (JHJ) is always diffraction-limited no matter how narrow the support of / is and 
how close the source is to the TRM. This means that in the absence of scattering even when 
the evanescent waves are recorded, phase-conjugated and back-propagated, the resulting field 
is still smeared below the scale of wavelength. In general, if |3?[G]| ^> l^fGJI in the near field 
as is the case for the free space, the latter does not contribute to the leading order effect of 
the near field TR and hence that TR, which affects only the phase, has a secondary effect 
on the subwavelength resolution of the near field TR. 

Since we are concerned here with the asymptotic behavior of near field TR as the point 
source approaches the TRM, we relax the definition of FWHM. Instead of "full width at half 
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maximum" in the literal sense we also refer to the (asymptotic) resolution as the asymptotic 
size of a focal spot as the maximum width at the same (or the minimum width at a slower) 
asymptotic behavior of the point-spread function. For example, in this metric, the focal spot 
size for (22) is A/2 while the truncated version of (23): 

cos (k\r — r'|) 



(25) 



mm 



P- 



47r|r — r' 



5 < 1 



has an asymptotic spot size 0(5). 

Next we shall analyze the near-field focusing property of the various TR schemes. We 
assume the near field asymptotic 

(26) |G(r,r')| ~ | 



r'l < A. 



For |r — r'| <C A (26) is the generally correct asymptotic as the near field singularity is 
determined solely by the Laplacian in the Helmholtz equation [19]. In addition to (26) we 
further assume 



(27) 



VG(r,r') 



r — r 



»/|3 : 



r'l < A. 



4.1. Transverse resolution of monopole TR. Consider the case of a point monopole 
source located at r in close proximity of A. Let r' Q be the closest point on A to r and 



r o 



5<1. Using (26) we have the estimate 

K M (r ,r ) = O(-klogS) 



since the singularity in the integrand of (10) is 

k 



(28) 



ro||r 



r - r -C A, 



fnl < A. 



More explicitly, the integral (10) can be simplified by replacing the surface integral over A 
by the tangent plane at ig. Without loss of generality, we may assume the tangent plane to 
be the x — y plane and r = (0, 0, 5). The leading order transverse profile of Km around r 
in x is given by 



(29) 



(x' 2 + y' 2 + 5 2 Y 1/2 {{x - x'f + y' 2 + 5 2 y 1/2 dx 'dy' ', 5 < 1 



where the integration is restricted to a bounded area, say the circular disk of radius A, around 
the origin to avoid logarithmic divergence at large x',y'. 

Let x = 5 a , < a < 1. Then (29) has the same asymptotic as 

(30) k J (x' 2 + y ,2 )- 1/2 ((5 a - x') 2 + y'^-^dx'dy', 5 < 1 

where the integration is over the punctured disk with holes of radius 5 centered at (0, 0) and 
(5 a ,0) since the two holes have a vanishing contribution to (29). The change of variables 
from x',y' to x" = x'/5 a ,y" = y'/S a renders (30) into 

(31) k [ (x" 2 + y" 2 )- 1/2 ((l - x") 2 + y" 2 Y l/2 dx"dy" 



Figure 2. (Left) 3d transverse profiles with FWHM = 
1.37(K M ),0.03(K D ),0.07(K MD ,K DM ); (Right) 3d longitudinal profiles 
with FWHM = 0.3565(ir M ),0.00421(ir D ),0.00947(if D M,^MD)- Note that the 
transverse profiles of the Kmd and K DM coincide. 

integrated over the punctured disk of radius \5~ a with shrinking holes of radius 5 1 ~ a centered 
at (0,0) and (1,0). Hence 

r , / \ ( —ka log 5, a G (0, 1) r 
(32) ^ M (r,r )~ 8 8«1 



where the logarithmic divergence for a > is due to the integration (31) on the scale of 
6~ a . 



When 5 is so small that log 6 is the dominant behavior of Km near the source point, (32) 
says that at x — 5 1 ^ 2 , Km is roughly at half the maximum and therefore the transverse 
resolution is proportional to 5 1 ^ 2 . However, this scaling behavior requires extremely small 
5 to manifest itself. Outside this regime, the focal spot size would decrease slowly as 5 
decreases (see the comments on Figure [2]). 

4.2. Transverse resolution of dipole TR. For r' = r' Q and r = r , r — r' is orthogonal to 
A and thus |VG| ~ \dG/dn\. Therefore the leading order transverse profile for Kn(r,r ) is 
given approximately by 

(<*q.) I / ^ x ^ dr'dn' 

1 ' k J {x> 2 + y> 2 + 8 2 f/ 2 {{x - x 1 ) 2 + y' 2 + 5 2 f' 2 y ' 



The integrand in (33) has a different singularity than that in (29). Indeed, for x = 5 a , the 



dominant contribution to (33) comes from the disks of radius 5 around (0,0) and (<5°,0) 
which is easily estimated to be 

K D (r,r ) ~ 5 1 - 3 , ae[0,l}. 

By contrast, K-£>(ro, ro) ~ 5~ 2 . Hence the focal spot size in this case is smaller than 5 a , Va G 
[0, 1) asymptotically. We conclude that the focal spot size for the dipole TR is essentially 
0(5) and that the dipole TR can in principle break the diffraction limit to achieve an 
arbitrarily fine resolution by bring the TRM sufficiently close to the source. 



Figure 3. The transverse resolution of various TR schemes in three (left) 
and two (right) dimensions. The linear behavior of R MD in the right plot 
demonstrates the square-root behavior of -Rmd in two dimensions. 



4.3. Transverse resolution of mixed mode TR. By the same analysis as before, the 



integrals (12) and (13) can be approximated by, respectively, 

5 



(34) 
(35) 



1 



(x 12 + y' 2 + 5 2 ) 1 / 2 ((x - x') 2 + y' 2 + <5 2 ) 3 / 2 
5 1 



dx'dy' 



x 



(x' 2 + y' 2 + 5 2 f/ 2 {{x - x') 2 + y' 2 + 5 2 ) 1 / 



Yj^dx dy . 



Set x = 5 a ,0 < a < 1. In the case of (34), the dominant contribution (~ 5 a ) comes from 



integration around (5 a , 0) for a < 1 while, in the case of (35), the dominant contribution (~ 
5~ a ) comes from integration around (0,0). In both cases i^MD(i"o, r o) ~ -^dm(i"o, r o) ~ 
Hence the focal spot size is smaller than 6 a , Va G [0, 1) asymptotically. 

We conclude that the focal spot size for the mixed mode TR is essentially 0(5) and that 
the mixed mode TR can achieve an arbitrarily fine resolution as the TRM approaches the 
point source. Moreover, the fact that the maximum values of the point-spread functions are 



inversely proportional to the focal spot size is consistent with the cut-off version of (23) in 
the following sense. Truncate the Green function (26) at the level 5 _1 as in (25). Then the 



asymptotic spot size of the resulting function is 0(5). In other words, the mixed mode TR 
of near field recovers qualitatively the asymptotic behavior of the Green function near the 
source point. 

When the superresolution terms from /Tdm and K MD combine in K PB , they cancel each 
other resulting in diffraction-limited focal spot. 

In Figure [2] the normalized transverse and longitudinal profiles are shown and their res- 
olutions calculated for various TR schemes. The profiles are normalized to be unity at the 
source location. In the simulation, we use A = 2tt, A = 2A and r = (0.01,0,0). Before 
normalization, the Km curve has the maximum less than 0.5 and hence is far from the 



square-root asymptotic regime discussed in the paragraph following (32). 
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4.4. Two-dimensional case. It is instructive to compare the scaling behavior of resolution 
in two and three dimensions. In two dimensions, (33) becomes 



(36) \ [ -J-^ x - * ■ dx'. 



k J x' 2 + 5 2 (x - x') 2 + 5 2 
The same calculation shows that 

K D (r,r )~5 1 - 2a 1 a G [0, 1] 
and K D (r ,r ) ~ The i^D-resolution asymptotic is given by 5 a ,a = 1. 



For the mixed mode, (34) and (35) become 



(37) lj i0g ^ + sr^x {x J )2+s ^ 

(38) -U x log ((x - x'f + 6*)-^dx' 



yielding 



K MD (v,r ),K DM (r,r ) ~ a (-log5) 2 , a E [0, 1] 



and i^ MD (r , r ), i^DM( r 0) r o) ~ ( — log 5) . Hence, the asymptotic resolution for the mixed 
mode TR is given by 5 a , a = 1/2. Empirically, the square-root regime for the mixed mode TR 
in two dimensions sets in much earlier than the monopole TR in three dimensions. Here lies 
the main difference between the two and three dimensional cases. Namely, in two dimensions 
the transverse resolution for the mixed mode TR is proportional to the square-root of the 
distance between the source and TRM. 

For the monopole TR in two dimensions, Km is uniformly bounded w.r.t. 5 and the 
resolution can not be deduced from the asymptotic analysis. It is likely though that the 
monopole TR resolution is qualitatively independent of 5. This is amply confirmed by our 
numerical simulations (Figure |3|. 

In Figure [3]the transverse resolutions, denoted by i?x>X = M, D, MD, DM, for the respec- 
tive TR schemes as a function of the distance between the source and TRM are plotted for 
the range [0.001,0.01]. The dipole resolution exhibits the linear behavior in two and 
three dimensions and the mixed mode resolution Rmd, Rdm exhibits the linear behavior in 
three dimensions and square-root behavior in two dimensions. Both the dipole and mixed 
mode resolutions are order of magnitude better than the monopole resolution. 

In Figure [4] the normalized transverse and longitudinal profiles in two dimensions are 
shown and FWHM calculated for various TR schemes with the same physical parameters as 
in Figure [2] The unnormalized transverse profile for the mixed mode TR has the maximum 
about 6. The dipole resolution is much better than the mixed mode resolution which in 
turn is much better than the monopole resolution in both the transverse and longitudinal 
directions. 

In summary, in the case of a monopole source the dipole TR produces the best resolution in 
both transverse and longitudinal directions, roughly in proportion to the distance between 
the source and TRM. The mixed mode TR can deliver nearly the same performance in 
three dimensions and roughly the square-root resolution in the transverse direction in two 
dimensions. Between the two mixed mode TR schemes, Kmd is still better than i^DM i n the 
longitudinal direction. A simple explanation for the superior resolution of TRM involving 
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Figure 4. (Left) 2d transverse profiles with FWHM= 
1.87(K M ),0M(K B ),0A7(K DM ,K MD ); (Right) 2d longitudinal profiles: 
FWHM= 1(K M ),0.009{K B ),0.U8{K MB ),0.22(K DM ). Note that the trans- 
verse profiles for the K MD and K DM coincide. 



the dipole field is that a dipole source has a stronger singularity than a monopole source (cf. 



(g7p versus g§) 

5. Approximate TR is biased 

We have already seen that the perfect TR described by ^ with G = G is diffraction- 
limited and contains no evanescent component of the initial outgoing field (see the remark 



following (24)). 

For the half space (domain (iii) of Figure [I]), the result is similar: Suppose the initial wave 



field is perfectly phase-conjugated on the TRM (i.e. the boundary condition (21) is satisfied 
exactly by, say, the scheme Kub)- Then in the interior of the half space the homogenous 
components of the initial wave are perfectly phase-conjugated but the evanescent components 
are exponentially damped by the factor 

e -2*v*?+*i-i, sj + sj > i 

depending on the distance z from the TRM to the field point (z, x) (Appendix [Cj) . 



On the other hand, for the paraxial wave (45), perfect phase conjugation on one plane 



implies perfect phase conjugation in the half space due to the fact that the paraxial wave 
field is uniquely determined by the direction of propagation and the initial condition on 
one plane (Appendix [C]) . As a consequence, vm and K M are perfect TR schemes in their 
respective contexts. In particular, Km can focus perfectly back on the location of the point 
source in the paraxial approximation. 

The perfect focusing of time reversed paraxial wave is, of course, unphysical because the 
evanescent waves are not accounted for. The other aspect that is more impractical than 
unphysical with the perfect conjugation of the homogeneous component of the spherical 
wave, is that an infinite aperture (A = M 2 ) is required for such a result. 

In this section, we show that when D is a half space and A is bounded, the maximum 
point of the TR focal spot does not correspond to the exact source location r in general. 
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In contrast, when A = IR 2 in the paraxial regime, since Km is a perfect TR scheme for the 
paraxial wave, the focal spot is always centered. 

This bias effect tends to be much more pronounced in the longitudinal direction. As the 
aperture of TRM increases, the transverse offset (bias) decreases while the longitudinal offset 
may persist if there is insufficient scattering in the medium, cf. Section |6j This is why we 
consider only the transverse resolution in the previous section. 

Let us recast the problem in more general terms. Let <p(r, r') be the phase of G(r, r') and 
define 

r 2 

(39) P M (r) = 



/ G*{r', r )e^ (r ' r,) ^(r, r')da{r') 

J A 



where W(r, r') is a non-negative function. The monopole TR kernel Km corresponds to 
W{r,r') = \G{r,r')\. 

Direct calculation of VPm(i"o) yields 

VPmM 

|G*(r\ r„)||G(r", r )|W(r , v')W{r Q , r")(iV0(r o , r ') - ^ V 0(r o , v"))da{v')da{v") 

AxA 

+ [ |G*(r',r )||G(r",r )|V(^(r,r')^(r,r")) r=ro ^(r')^(r"). 

J AxA 

The first term vanishes due to anti-symmetry with respect to r', r". 

Consider the case of a planar TRM on the plane z = with the paraxial wave propagating 
in the ^-direction (Appendix A). Suppose the medium inhomogeneities act like a phase object 
affecting only the phase of the free propagator. Namely, G(z, x, z',x') takes the form 



■ k |x — x 



'|2 



(40) G(z, x, z', x') = iHz-z^j^^rjjeiz^z'^) ^ £ r2 _ 

v ' 4:ir(z - z') 



Note that for (40) the expression (39) is divergent and thus invalid at r = ro if A = K . 

For a finite aperture, W(r, r') = (4ti\z — z'l)" 1 and thus V (W(r, r')W(r, r")) r=ro points 
in the negative z direction and consequently ro tends to be farther away from A than the 
maximum point of Pu- In other words, the focal spot of the monopole TR tends to be closer 
to A than the source location. Even in the case of a closed cavity the offset (bias) is reduced 
but not completely eliminated except for special cases. 

Now let us change the perspective and think of TR as a means of imaging a point source 
in a perfectly known medium. In this context a few widely used W(r, r')'s are: 

W(r,r') = IG^rOr 1 (the inverse filter), 

W(r,r') = W(r')/J J \W{r")\ 2 da(r"), W G L 2 {A,da) (the phase processor), 

W(t,t') = \G(r, r ')\/ J J \G(r,r")\ 2 da(r") (the conventional matched field processor). 

For the inverse filter, V(W(r, r')W(r, r")) r=ro and thus VPm(i"o) tend to point away from 
the TR mirror A and hence the source location tends to be closer to A than the maximum 
point of Pm- Again, when A = 3D or in the far field regime this offset can be reduced. 
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Direct substitution shows that VPm( f o) — for the two latter choices of W. Moreover, 
r is the maximum point of Pm for both choices of W. In the case of the W-weighted phase 
processor, we clearly have 



Pu(r) < 



\G*{r',r )\W{r')da(r') = P M (r ). 



In the case of the conventional matched field processing, it is well known [21 EH] that 



(41) V(r , r') = G(r , v')/J jf |G(r , r') \ 2 da(r>) 

is the solution to the optimization problem: Maximize the quantity 

2 

G*(r',r )V(r')da(r') 



over V subject to the constraint 
(42) 



J \V\ 2 (r')da(r') = 1. 



Therefore 



P M (r)= / G*(r',r )V(r,v')da(r') 

J A 



is maximized at r = ro with the conventional matched field processing. 

In summary, as an imaging method both TR and the inverse filter are generally biased; the 
former tends to underestimate the range of the source while the latter tends to overestimate 
it. The conventional matched field processor is unbiased and produces the optimal signal- 
to-noise ratio. The phase processor, which uses only the phase of the Green function, is a 
simple alternative to the conventional matched field processor. 



Moreover, in the regime described by (40), \G(r, r')| is a function of the distance between 
the source point and the TRM. In such a case the conventional matched field processor 
coincides with the phase processor with W(r') = Ia/IA] 1 ^ 2 independent of r'. 

Similar conclusions can be drawn for the pure dipole case with 



(43) 



^o(r) 



dG*(r',r ) 
dn 



W(r,r')dcr(r') 



where if>(r, r') is the phase of dG(r, r')/dn. Pd can be analyzed as above with G replaced by 
its normal derivative. 



For Kmd, K-£,m, consider the paraxial regime for which (19) is valid. Each term on the 



right hand side of ( 19 ) can be analyzed as before. The result is that the mixed mode TR's 



are biased, especially in the longitudinal direction, just like the monopole and dipole TR's. 



6. Randomness eliminates bias 

In this section, we explain how randomness can make the TR focal spot centered at the 
location of the source point. 

Consider typical profiles of the Green function in a multiple scattering medium as de- 
picted in Figure |5j A main feature of the intensity profiles is the rapid fluctuations in the 
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1000 1500 2000 2500 3000 3500 4000 4500 5000 



Figure 5. The transverse (left), at x — 5000, and longitudinal (right), at 
y = 2500, profiles of the intensity of the 2 — d Green function at wavelength 70 
with 1000 particles randomly distributed in [2000,4000] x [0,5000] and source 
(not shown) at (—5000,2500) whose n 2 — 1 is 70. (Reproduced from [TB] ) 



transverse and longitudinal directions. These fluctuations upon differentiation and averaging 
(integration over A x A) as in the expression 



(44) 



VP M (r ) = \j V\G{v^)G{vy)\ 2 r=rQ da{v')da{r'') 



essentially cancel one another. 

More precisely, the degrees of scattering and randomness can be measured by the spatial 
spread of the propagation channel. If the spatial spread increases to infinity, then the TRM 
acquires an effective aperture that is proportional to the spatial spread jU [T3] . Then, as per 
discussion at the beginning of Section |5j the accuracy of approximate TR's improve and the 
focal spot becomes more centered and sharper. 

Moreover, counterintuitively, the TR image can avoid statistical fluctuations if the medium 
is sufficiently random. Let us recall some analytical results from the literature. 

Consider TR in a random medium occupying a half space with TRM on the planar bound- 
ary. When the random channel is the Rayleigh fading regime (i.e. obeying zero-mean Gauss- 



ian statistics) with a divergent spatial spread the monopole TR (10) with TRM's aperture 
much larger than the coherence length i c will focus on the source point with a spot size 
~ i c [T3J E]. The conventional coherence theory suggests that the coherence length is of 
the order of wavelength [23]. However, spatial correlations of much smaller extent has been 
recently measured in the near field of random media [I]. This effect in conjunction with the 
preceding theory may be the key to understanding the subwavelength focusing observed in 
the time reversal experiment reported in |21j . 

The result of stability and focusing is a special case of what was originally established for 
the setting of multiple frequencies and multiple source points [T3| III]. With straightforward 
modification of the arguments in JT3J [H] the same result can be extended to any TR of the 



types ( 10 )-( 13 ) in the paraxial regime, the only difference being in the asymptotic shape 
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of the focal spot. Experiments show consistently similar results on stability and focusing 
properties [171 1211 127] . 

In conclusion, random fluctuations in the intensity of Green's function is a necessary 
condition for the desired focusing property of time reversal with a finite aperture. 



7. Conclusion 

We have reexamined various time reversal schemes that can be employed from the bound- 
ary of a domain; they involve monopole and/or dipole fields. 

We have explored their relationships in the paraxial regime. For the boundary behavior, 
we show that for paraxial wave only the monopole TR scheme satisfies the exact boundary 
condition while for the spherical wave only the MD-mode TR scheme satisfies the exact 
boundary condition. We have also shown that for the paraxial wave the standard monopole 
TRM produces the perfect result in the entire domain but for the near field the monopole 
TRM can hardly achieve subwavelength focusing. On the other hand, the dipole TRM is 
approximate in the paraxial regime but is capable of producing a focal spot size linearly 
proportional to the distance from TRM to the source location, thus focusing on a subwave- 
length scale when the TRM is sufficiently close to the source point. The mixed mode TRM 
has the similar (linear) behavior in three dimensions but the square-root asymptotic in two 
dimensions. The monopole TRM is, if possible at all, ineffective to focus below wavelength. 

We have seen how two mixed mode effects combine in the so called perfect TR scheme 
for a closed cavity and cancel their respective subwavelength focusing property. The perfect 
TRM is "perfect" only for a source-free initial field. 

We have examined an often neglected effect of biased focusing associated with TRM of 
a finite aperture in a weak- or non-scattering medium. This effect pertains to all five TR 
schemes discussed in the paper. In a multiply scattering medium, the bias is greatly reduced 
and the focal spot becomes centered at the source location. The removal of bias by random 
scattering is attributed to the enlarged effective aperture and the random fluctuations in the 
intensity of the Green function. 



Appendix A. Mixed mode TR for paraxial wave 

Let z G K be the longitudinal coordinate and x G R 2 the transverse coordinates. Let 
Gp(z, x, z 1 , x') be the Green functions for the paraxial wave equation 

(45) ± i^-V + ^-A T ^ + k -[n\z, x) - 1)$ = 

az 2k 2 

where At is the transverse Laplacian. The plus sign represents positive z propagating wave 
and minus sign negative z propagating wave. 

The full Green function G(z, x, z', x') in the paraxial regime can be expressed as 

(46) G(z, x, z 1 , x') = — — Gp(z, x, z',x!), z > z 

(47) G(2,x,£',x') = 16 2k G^fox^x'), z<z' 

Let D = {z > 0}, the plane z = be occupied by the TRM and r = (-2 ,x ) be the 
location of the source. Assume that the TRM is away from the medium inhomogeneities i.e. 
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n(0,x) = 1. Then K MD can be written as 



(48) K UD {r,r ) 



+ 



4k 

e ik{z-z ) 



G P *(0, x', Zo, x )Gp (0, x', z, x)dx' 



4k 2 



f 9 

/ G P * (0, x, zo, x )—G P (z , x, z, x) | 2 / =0 rfx'. 



By the assumption n(0, x) = 1 and (45), we have 



(49) 



K MT )(r,r ) 



+ 



j^ e ik(z-z ) 

4k 

i e ik(z-zo) 



Ik 3 



G P * (0, x', zq, x )G P (0, x', z, x)dx' 
G P *(0, x', z , xq) A' T Gp (0, x', z, x)dx' 



— i/c y G*(0, x', zo, x )G(0, x', z, x)dx' 

+ ^ / V^(0,x',2„,x ). V^ r G(0,x , ,z,x)dx' 



after integration by parts. The result for Kdm can be likewise derived. 



Appendix B. Boundary values 



B.l. Paraxial wave. Let us first focus on the boundary values of (14)-(17). Let w(z,x) be 



a paraxial wave propagating in the negative z direction which can be written as 



e- ife ^(z,x) 



where \I/ is the solution of eq. (45) with the minus sign. 



In the paraxial regime with n(0, x) = 1, we have from (46) and (45) that 



«m(0,x) 
ud(0,x) 
w M d(0,x) 
v D m(0,x) 



u*(0,x')G£(0,x,0,x')<£ic / 
y A T u*(0,x')G+(0,x,0,x')dx' - ^ y A^«*(0,x / )G+(0,x,0,xOt2x' 
^ y n*(0,x , )G+(0,x,0,x , )dx' + ^ J A T u*(0, x')G+(0, x, 0, x')(2x' 
^ M*(0,x , )G+(0,x,0,x / )rfx' + ^y At^CxOG+^x^x'^x' 
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where u is a paraxial wave propagating in the negative z direction. Setting Gp (0, x, 0, x') 
5(x — x'), we obtain the following boundary conditions 



um(0,x) 


= u*(0,x) 


f D (0,x) 


= --^A T u*(0,x')-^A^* 


w M d(0,x) 


= ^*(0,x) + i ^A T «*(0 ) x) 


w D m(0,x) 


= i«*(0,x) + ^ 5 A T ti*(0 J x) 


D PB (0,x) 


= u*(0,x) + ^A T u*(0,x). 



The kernel ii'x, X = M, D, MD, DM, PB satisfy the similar boundary properties. In other 
words, only the monopole TR gives the correct boundary values in the case of paraxial wave 
propagating in the half space. 

B.2. Spherical wave. Consider the half space D = {z > 0} with A = {z = 0}. Assume 
G = Gq, the free space Green function. 

Differentiating the Weyl representation for Gq with respect to z', we find that 



(50) 



d_ 

dz 



7 G (z,x, z',x' 



\z=z'=0 



-5(x 



From (50) it is immediately clear that 

1 



1 



w M d(0,x) = -m*(0,x), K ME) (0,x) = -G*(0,x,2 ,x ). 

However, neither vbm nor Kdm satisfies the boundary condition as we show now. 
In the neighborhood of A, the wave field u in D can be represented in terms of the plane 
wave spectrum: 

(5lXx,y,z) = J ds 1 ds 2 e ik{siX+S2y) e lkS3Z A(s 1 ,s 2 ) + J ' d Sl ds 2 e ik{siX+S2y) l e~ ikssZ B( Sl , s 2 ) 
where S3 is given in (24). Consider for simplicity one-way wave with, say, 5 = 0. Then the 

d Sl ds 2 e ik{siX+S2y) e- lkssZ A*(- Sl ,-s 2 ) 
ds 1 ds 2 e ik ^ slx+S2 ^e ikssZ A*{-s h -s 2 ) 



phase-conjugated field has the form 

u*(x,y,z) = 



sf+s|<l 



+ 



i+si>l 



and thus 
d_ 

dz 



u*(x,y,0) 



k / d Sl ds 2 e 4fc(si * +S2 ^ 3 ^*(-si, -8 2 ) [-/ {S?+S 2< 1} + / {s?+s | >1} 



where Is is the indicator function of the set S. By Plancherel's identity we obtain 

-I{si+s 2 2 <i} + h4+4>i} 



w DM (0,x) = 2vr y ds 1 ds 2 e ik{siX+S2y) A*(-s 1 ,-s 2 ) 
^ 2vra*(0,x). 



Namely, no constant multiple of vdm satisfies the boundary condition. 
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Likewise, one can show that none of the other TR schemes than v md, ^md satisfies the 
boundary condition, with some adjustable constant factor. 



Appendix C. One-plane TRM in a half space 



Consider the half space D = {z > 0} with A = {z = 0} and point source located at 
(zo,x ), zq > 0, Figure [6j We assume that the medium inhomogeneities n 2 — 1 is away from 
the TRM so that the gap D~ is a non-empty infinite slab. Figure [6] 

First we observe that for the paraxial wave, perfect phase conjugation on z = implies 
perfect phase conjugation throughout the domain D. This can be shown as follows. 

Let u(z,x) = e~ lkz ^~ (z, x) be the initial wave propagating in the negative z direction 



where ^ satisfies (45) with the minus sign. If m(0, x) is perfectly conjugated by, say, either 



the scheme described by i>m or t> PB , then the resulting field, denoted by v(z,x), is given by 



u(z,x) = e 4fe ^+(z,x) 



where satisfies (45) with the plus sign and the initial condition 

* + (0,x) = tf-*(0,x). 



Now since \P *(z,x) satisfies (45) with the plus sign and the initial condition \l/ + (0,x), we 
have from the uniqueness theorem for the initial value problem of (45) that 



and consequently 



V\Z, X 



Z,X) 



u 



z>0. 



The same argument and conclusion apply to the paraxial wave with a point source. 

The case with the spherical wave is similar but more subtle since the evanescent wave is 
involved. It can be shown that the Green function in D~ n {z < z } can be expressed as 
superposition of plane waves propagating in the negative z direction 



(52) 



G(r,r ) 



ds\ds2e 



ik(siX+S2y) —ikssz 



B(si,s 2 ), r E D f]{z<z } 



in analogy to (51 ). 



Let W be a source-free field propagating across z = into D and be represented as 
(53) W(t) = J d Sl ds 2 e lk{siX+S2y) e ikS3Z A( Sl ,S2), r £ D~. 



Note the sign in the exponent associated with z here is different from that in (52) since 



they propagate in the opposite directions. Suppose W is the phase conjugate field of G at 
z = 0, i.e. W*(x,y,0) = G(x,y,Q,r ),Vx,y. For z = 0, equating (53) with (52) conjugated 
we obtain the condition 



(54) 



B* 



•si, -s 2 ). 



This result is first observed in [22] . 

Let Wh and We be, respectively, the homogeneous and the evanescent components of W, 
corresponding to, respectively, the integration (53) restricted to s\ + si < 1 and s\ + > 1. 
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Figure 6. The dark area represents the scattering domain S. The asterisk, 
*, represents the source point r which may be anywhere in D. 



The homogeneous and the evanescent components Gh,Gi of (52) are defined analogously 
Then (54) lead to, respectively, 

(55) W H {r) = G„(r), r e Br 

(56) W K r) = JM^^V^BM, reir. 

In other words, the homogenous components of the initial wave is perfectly phase-conjugated 
but the evanescent components are exponentially damped by the factor 

e -2kZy/s\+sl-l 

depending on the distance from the TRM to the field point. 
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